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I .  INTRODUCTION 


A yawsonde  is  a device  developed  to  give  a continuous  record  of 
the  angle  between  a shell's  axis  of  rotation  and  the  line  of  sight  to 
the  sun1-6.  It  consists  of  two  light  sensors  located  behind  two 
inclined  slits  on  the  shell  body  and  a telemetry  unit.  The  instant 
each  sensor  sees  the  sun,  a discrete  pulse  is  transmitted  to  a ground 
station  and  recorded  versus  time.  The  resulting  pulse  train  can  then 
be  analyzed  to  yield  sun  angles  to  an  accuracy  of  0.1°  and  the  spin 
rate.  The  usual  reduction,  however,  assumes  the  sun  angle  to  be  con- 
stant for  a single  revolution  of  the  shell.  If  the  shell  is  perform- 
ing large  high-frequency  angular  motion,  this  assumption  can  be  invalid 
and  erroneous  data  can  be  produced.  It  is  the  purpose  of  this  report 
to  develop  a valid  data  analysis  system  to  handle  yawsonde  data  without 
assuming  the  constancy  of  the  sun  angle  during  a single  revolution  of 
the  shell. 


II.  CONSTANT  SUN  ANGLE  ANALYSIS 

We  will  make  use  of  a solar  fixed-plane  coordinate  system  which 
can  be  described  in  terms  of  three  unit  vectors,  ej  is  along  the 

missile's  axis,  e3S  is  in  the  plane  determined  by  the  sun  and  the 


1.  I.  Amery , H.  Henning , K.  Laurie,  and  E.  Wlatnig,  "A  Telemetry 
System  for  the  Measurement  of  the  Yaw  of  a Projectile  Throughout 
the  Major  Part  of  Its  Trajectory  (U) ,"  RABDE  Report  1/65 , March 
1965.  AD  359250.  (Confidential  Report) 

2.  W.  R.  Haseltine,  "Yawing  Motion  of  5.0"  MK  41  Projectile  Studied 
by  Means  of  Yaw  Sondes, " Haval  Weapons  Center  Report  NWC  TP  4779, 
August  1969.  AD  862065. 

3.  W.  H.  Mermagen,  "Measurements  of  the  Dynamical  Behavior  of  Pro- 
jectiles over  Long  Flight  Paths,"  Journal  of  Spacecraft  and 
Rocke ts,  Vol.  8,  April  1971,  pp.  380-385.  (See  also  BRL  MR  2079. 
AD  717002) 

4.  W.  H.  Clay , "A  Precision  Yawsonde  Calibration  Technique,  " BRL  MR 
2263,  January  1973.  AD  758158. 

5.  R.  H.  Whyte  and  W.  H.  Mermagen , "A  Method  for  Obtaining  Aerody- 
namic Coefficients  from  Yawsonde  and  Radar  Data, " Journal  of 
Spacecraft  and  Rockets , Vol.  10,  June  1973,  pp.  384-388.  (See 
also  BRL  MR  2280.  AD  759482) 

6.  W.  H.  Mermagen  and  W.  H.  Clay,  "The  Design  of  a Second  Generation 
Yawsonde,"  BRL  MR  2368,  April  1974.  AD  780064. 
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missile's  axis  and  points  toward  the  sun's  side  of  the  missile's  axis, 
while  ?2s  is  t3s  X lx. 

Each  sun  sensor  unit  emits  a pulse  at  the  instant  that  the  sun  is 

in  the  plane  of  its  inclined  slit.  Let  y be  the  inclination  of  the  first 
1 a 

slit  with  respect  to  the  missile's  axis  and  let  <J>a  be  the  missile's 

roll  angle  with  respect  to  the  vector  £2s  when  the  line-of-sight  from 

the  first  sensor  to  the  sun  lies  in  the  e.e  plane.  Then  the  vector 

1 38  r 

normal  to  the  plane  of  the  first  sun  sensor  is 

N = sin  y + [ t2s  cos  + ^3s  sin  (<MJ  ] cos  y (1) 

d a a a a 

where  <p  is  the  roll  angle  of  the  missile.  The  sun  line-of-sight  can 
be  specified  by  the  unit  vector 

s = ex  cos  a + £3s  sin  a (2) 

where  a is  the  angle  between  the  sun's  rays  and  the  missile's  axis. 

When  the  sun  is  in  the  plane  of  the  first  sensor 

ft  . s = 0 (3) 

a 

or 

sin  (<J>-<|>  ) = - tan  y ctn  a (4) 

a a 

A similar  relation  applies  to  the  second  sun  sensor: 

sin  O-^)  = - tan  ctn  a (5) 

We  now  consider  a triple  of  time  pulses  consisting  of  two  sight- 
ings of  the  sun  by  the  first  sensor  at  t = tx,t3  and  one  by  the  second 

sensor  at  t = t2-  The  pairs  of  <|k,  ck's  for  these  three  times  satisfy 

the  relations: 


sin  ) = 

- tan  y 

ctn  a. 

i = 1,3 

(6) 

l a 

a 

l 

sin  02-<f>b)  = 

- tan  yb 

ctn  a2 

(7) 

We  now  assume  that  0^  = 02  = 03  = 0 and  that  | is  a constant.  Then 

sin  = sin  (<f>3-<J>  ) (8) 

a a 
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so  that* 


or 


Now 


$3  ~ $1  - 2tt 


(9) 


<j>  = 


2tt 

t3_tl 


CIO) 


(f>2  - ((>1  = 4>  Ct2”ti) 


= 2tt  Ct2-ti)  (t3-t1)  1 


(ID 


Let 


<S  = <j>2  - <h  _ (<£b  - <f>a) 


(12) 


Then 


sin  (<p 2 - <|>b)  = sin  (<i> x - <j>a  + <5) 

= sin  (<j>i  - <j>  ) cos  6 + sin  6 cos  (4)!  - <j>  ) (13) 

3.  3- 

Equations  (6-7)  can  now  be  substituted  in  Equation  (13)  and  the  result 
simplified  to  yield 

tan  a = [tan2  y - 2 tan  y tan  y,  cos  6 + tan2  y,  ] 2 (sin  6)-1 
a a b D 

(14) 

This  equation  gives  the  sun  angle  in  terms  of  the  three  times, (t^,  t^, 

t3)  and  three  parameters  of  the  sensors  (ya»  y^,  4>b-4>a)  and  Is  valid 

for  idealized  sun  sensors.  Actual  sun  sensors  will  not  exactly  satisfy 
Equations  (4-5)  but  roll  angles  for  each  sun  sensor  can  be  measured  as 
functions  of  the  sun  angle4. 

♦ . = f (a.)  + (i-1)  it  i = 1,  3 (15) 

l a i 

<P2  = f (a2)  + 2t ;l  (16) 


* Since  the  sun  sensor  plane  is  really  a half  plane,  the  root 
<j>3-<j>i  = Tr  is  not  appropriate. 
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where  l is  an  integer  selected  so  that 


Once  again  we  assume  that  03=02  = 03  = 0 and  that  ^ is  a constant. 
Equations  (10-11)  are  still  valid  and 


<t>2  - <h  = fbO)  - faC<0  + 2tt j = g(a) 


(17) 


or 


a = g 1 (<|>2 -<!>!) 

where  1 s given  by  Equation  (11). 


(18) 


III.  VARYING  SUN  ANGLE  ANALYSIS 

We  assume  that  the  sun  angle  is  a linear  function  of  time  and  that 
the  roll  angle  is  a quadratic  function  of  time  and  hence  we  need  five 
successive  sightings  of  the  sun  to  determine  the  five  parameters: 

o±  = o 3 (ti-t3)  + a3  i = 1,2, 3, 4,5  (19) 

♦ = <J>3  (t^-t2)2/2  + $3  (t^-t2)  + <j>3  (20) 

If  three  of  the  sightings  are  by  the  first  sensor  and  two  by  the  second, 
then 


= f (a.)  + (i-l)ir  i = 1,3,5  (21) 

1 d 1 

4> . = f (a.)  + 2irZ  + (i-2)  tt  i = 2,4  (22) 

l b l 

The  five  parameters  (o3,  a 3,  4>3,  $3>  <j>3)  ,ean  be  obtained  from  the  five 

times  (t*,  t2,  t3,  t4,  t5)  by  substituting  Equations  (19-20)  in  Equations 

(21-22)  and  solving  the  resulting  nonlinear  algebraic  equations  by  an 
iterative  differential  correction  technique7.  The  first  set  of  values 
of  (a3,  03),  i.e.  (o31,  033)  can  be  obtained  by  using  the  constant-sun- 

angle  solutions  for  (tj,  t2,  t3)  and  (t3,  t^,  t5) . Since  these 


7 . P.  R.  Bevingtorij  Data  Reduction  and  Error  Analysis  for  the 
Fhysical  Sciences , MeGraw-Htll,  New  York , 1969. 
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approximately  correspond  to  the  sun  angles  for  t2  and  t4,  we  will 
denote  them  by  02,04. 

cr31  = (a2  + 01+3/2  (23) 

031  = (a4  - a2)/(t4  - t2)  (24) 


Equations  (19-20)  for  the  (j-1)  iteration  can  be  written  in  the  form: 


*3j  ^i  ' t3)2/2  + ^3 j (^  * *3)  + ^3 j 


= f„  (o. . 

K v ij 


1-F  +( 


8 fK 


3o3 


a3i-i-) 


3 j 33 


3 fK 


3a  1 


(a3j  " °3j-^  + miT 


(25) 


where 

K = a,  m = i-l  for  i = 1,3,5 
K = b,  m = 2(1-1)  + i for  i = 2,4 

Equations  (25)  are  linear  in  the  unknowns  (0,.,  a , 4.,.,  4 .,  4 .)  and 

3j  3j  OJ  33  3 J 

can  be  easily  solved;  usually  only  three  or  four  iterations  are  required. 


A 

can  be 
by: 


second  version  of  a varying  sun  angle  analysis  using  five  times 
used  for  constant  <J> . In  this  case.  Equations  (19-20)  are  replaced 


a. 

1 


o3  (t.  - 1 3)  2/2  + O3  (t.  - t3)  + a 3 


(26) 


♦i  = ^3  Cti  - t3)  + 4 3 (27) 

Equations  similar  to  Equations  (25)  can  be  easily  derived  for  this 
variant . 


IV.  EFFECT  OF  MOVING  COORDINATE  SYSTEM 

— v • 

41  is  the  roll  angle  measured  with  respect  to  e2S-  4 is  a roll  rate 

measured  in  this  solar  fixed-plane  coordinate  system  and  is  not  necessarily 
equal  to  the  roll  rate,  p,  measured  in  an  inertia  coordinate  system  since 
the  solar  coordinate  system  can  itself  have  a non-zero  roll  rate. 
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For  our  inertia  system,  we  will  use  an  earth-fixed  coordinate 
system  Xj,  X2,  X3,  so  aligned  that  the  X3  axis  is  vertical  and  the 
trajectory  is  initially  in  the  X1-X3  plane.  Since  the  actual  trajectory 
can  be  reasonably  well  approximated  by  a planar  trajectory,  we  will 
approximate  the  velocity  vector  along  the  trajectory  by  (V  cos  0T,  0, 

V sin  0T) . 1 

In  the  appendix,  the  roll  component  (tols)  of  the  angular  velocity 

of  the  sun  fixed-plane  coordinates  is  computed  in  terms  of  the  hori- 
zontal and  vertical  components  of  the  complex  yaw  in  the  aeroballistic 
fixed-plane  coordinates®: 


where 


“is  * 0>i  ?H  + b2  Cy)  COS  aT 


+ ^3  ?U5U  + bi*  + bs  + bg  ?u5u  (28) 


i 


= sin  6 + i cos  6 sin  a 


a,  0 are  the  angles  of  attack  and  sideslip 

aT  is  the  angle  between  the  sun  vector 
1 and  the  tangent  to  the  trajectory; 

b. 's  are  defined  in  Table  1. 


TABLE  1.  COEFFICIENTS  IN  EQUATIONS  (28)  AND  (A17) 
bj  = (s3  cos  eT  - sl  sin  0T)  esc2  oT 

b2  = -s2  esc2  ctt 


= b2  [2  cos 

°T 

tan 

0T 

- (1 

+ 

cos2  aT) 

bll 

= b2  [2  cos 

°T 

tan 

0T 

+ (1 

+ 

cos2  oT) 

= cot2  oT  - 

bl 

cos 

oT 

tan 

0T 

- b22  (1 

+ cos2  aT) 

bg  = -cot2  - bj  cos  aT  tan  0T  + bj2  (1  + cos2  oT) 


8.  C.  H.  Murphy,  "Free  Flight  Motion  of  Symmetric  Missiles", 
BRL  Report  1216,  July  1963,  AD  442757. 
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The  variation  of  the  complex  yaw  can  be  described  by  the  usual 
epicycle 

U,  i <P2 

5 = K1  e + K2  e 

- 4>jo  + t (29) 

The  presence  of  derivatives  in  all  terms  of  Equations  (28)  suggests 
that  the  high  frequency  mode  is  the  more  important  mode  in  Equation 
(29).  If  we  neglect  the  low  frequency  mode  (K2=0) , Equation  (28) 
becomes 

“is  = t_bl  Sin  *1  + b2  C0S  *1  Ki  cos  °t 

+ [(b4  - b3)  sin  <j>1  cos  <j>  1 + b5  cos2  4>1  - bg  sin2  4^]  4^ 

(30) 

Since  4>  is  the  difference  between  the  projectile  spin  and  the 
coordinate  system  spin: 

41  s P - (31) 

it  will  be  a periodic  function  with  frequency  <J>1  except  when  = tt/2. 
For  this  special  case  where  the  sun  is  perpendicular  to  the  projectile, 
the  frequency  will  be  2^.  The  average  value  of  <p  is,  however,  quite 
close  to  the  projectile  spin: 

ilv'P*  (>>5  - V *1  Kl2/2  <32> 


V.  DISCUSSION 

In  Figure  1,  the  sun  angle  complement  (a  = tt/2  - a)  and  <p  as 
obtained  from  the  three-time-measurements,  constant-sun-angle  reduction 
are  plotted  for  an  eight-inch  projectile.  These  data  are  then  reduced 
by  the  f ive-time-measurements  analysis  of  this  report  and  plotted  in 
Figure  2.  The  sun  angles  are  changed  by  very  little  but  the  oscillations 
in  $ are  reduced  by  50%.  The  remainder  is  clearly  periodic  with 
frequency  4^.  Later  in  the  flight,  the  sun  angle  went  through  tt/2; 

Figure  3 shows  a and  4>  for  this  portion  of  the  trajectory.  We  note 
• • 
the  4>  oscillations  are  much  smaller,  with  frequency  2<P1. 
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It  is  interesting  to  note  that  the  need  for  the  varying  a data 

reduction  can  be  eliminated  by  a modified  yawsonde.  If  the  first  sun 

sensor  is  oriented  so  that  Y * 0,  that  sensor  is  insensitive  to  sun 

a 

angle  (see  Equation  (4)),  so  that 

f (a)  = 0 (33) 

cl 

This  sensor  directly  measures  41  . The  second  sensor,  then,  yields  a 
a2  at  t=t2  from  the  standard  three- time-measurements  procedure. 
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Figure  1A.  Sun  Angle  versus  Time  for  Eight-Inch  Shell  Computed  from  Sets  of  Three  Sun  Sights 
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24 


OF  SYMBOLS  (CONTINUED) 

A A 

+ i £v  , the  complex  yaw  in  the  aero- 
ballistic  fixed-plane  system 

sin  B,  the  horizontal  component  of  £ 

a a A 

cos  6 sin  a,  the  vertical  component  of  £ 


the  sun  angle:  the  angle  between  the  sun 
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APPENDIX  A.  DERIVATION  OF  EQUATION  (28) 

Our  earth-fixed  coordinate  system  X,,  X2,  X3  will  be  so  oriented 
that  the  X1  - X3  plane  is  the  vertical  plane  containing  the  initial 
velocity  vector.  If  we  make  the  very  reasonable  approximation  that 
the  velocity  vector  stays  in  this  plane,  it  can  be  written  in  the  form 

^ = V (cos  0T,  0,  sin  0 ).  (Al) 

The  unit  vector  along  the  missile's  axis  then  takes  the  form 

e1  = (n3,  n2,  n3)  = (cos  0T  + m>  n2,  sin  0T  + n3)  (A2) 


nl2  + n22  + n32  = 1 


(A3) 


where  the  n ^ 1 s become  zero  for  zero -amplitude  yawing  motion. 


We  will  be  using  two  coordinate  systems  which  use  two  different 
fixed  planes : 

1.  The  usual  aeroballistic  fixed  plane  which  contains  the 
missile's  axis  and  the  vertical  vector; 

2.  A sun  fixed  plane  which  contains  the  missile's  axis  and  the 
sun  vector  and  is  the  reference  plane  for  the  yaw  sonde's  4> . 

The  other  two  unit  vectors  for  the  aeroballistic  fixed-plane 
system  are 


(0,0,1)  x e1  (-n2,  cos  eT  + n3,  o) 


'2a  cos  0 


cos  0 


(A4) 


x [(0,0,1)  X e3]  (0,0,1)  - (sin  0T  + n3)  e1 


3a 


cos  0 


cos  0 


(A5) 


cos  0 = [cos2  ©T  + 2nj^  cos  0T  + n32  + n22] 2 


(A6) 
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Similarly,  the  other  two  units  vectors  for  the  sun  fixed-plane  system 
are: 


-*•  ■+ 

- s X ei 

e = — 

2S  sin  a 


(A7) 


X (s  X e 1 ) s - (cos  a)  e1 


' 3s 


sin  a 


sin  a 


(A8) 


where 

s = (si,  s 2 » S3)  is  the  unit  vector  pointing  to  the  sun, 

-y 

cos  a = s • e, , 

and  a is  in  the  first  quadrant  when  the  sun's  rays  illuminate  the 
missile's  nose  and  is  in  the  second  quadrant  when  the  sun's  rays 
illuminate  the  missile's  base. 

The  horizontal  and  vertical  components  of  the  complex  yaw  in  the 
aeroballistic  fixed-plane  coordinates  can  be  computed  from  their 
definitions : 


«H" 


e2a  ' V "n2  C0S  0T 
V = cos~e 


(A9) 


= 


e_  • V 
3a 


n1  sin  eT  cos  0T 


n3  (1  + sin2  0T  + n : 
cos  0j 


cos  0T 


nq  sin  0 ) 
L_  (A10) 


where 

£ = £ + i^  = sin  6 + i cos  6 sin  a 

H V 

and  6,  a are  the  angles  of  attack  and  sideslip  respectively.  A 
quadratic  approximation  for  m and  a linear  approximation  for  cos  0 
can  be  obtained  from  Equations  (A3,  A6,  A9-10). 


28 


nx  = -n3  tan  e 


(All) 


n22  cos2  eT  + n32 


2 cos 3 en 


cos  e 
cos  0 


l + 


n3  sin  0T 

COS2  Qrr, 


(A12) 


n2  = -^H  + tan  0T 


(A13) 


A A 

n3  = -Cv  cos  0T  + (Cv2“Ch2) (sin  0t)/2 


(A14) 


The  roll  component  of  the  angular  velocity  of  the  sun  fixed-plane 
coordinates  can  now  be  computed  in  terms  of  the  n^'s: 


is 


= -e 


2s 


e 


3S 


= -e2g. [(esc  a) 


s - (ctn  a)  * e - ctn  a ej 


= (ctn  a)  e2s  • e3 


= ctn  a esc  a [n1  (-s2  sin  0^  + n2  s3  " n3  s2^ 

+ n2  (s3  sin  0T  - s3  cos  0T  + sL  n3  - s3  n3) 

♦ n3  (s2  cos  0T  + s2  ni  " si  n2)  (A15) 


cos  a 


1 


= cos  aT  + sL  nx  + s2  n2  + s3  n3 


(A16) 


where 


cos  aT  = si  cos  ©T  + s3  sin  0T 


29 


We  can  now  use  Equations  (All- 13,  A16)  to  express  u)ls  as  a 
quadratic  expansion  in  and 


"is  * (bl  «H  * b2  EV>  C0S  °T 

• • • • 

* b 3 VH  * \ Mv  * b5  Vv  * b6  EHEV  (A17) 
where  the  b^'s  are  defined  in  Table  1. 

Equation  (A17)  is  precisely  Equation  (28)  of  the  text.  It  is 
interesting  to  note  that  for  the  special  case  of  the  sun  directly 
overhead,  the  two  fixed-plane  coordinate  systems  are  the  same.  For 
this  case, 


s 


1 


= 


0 


flly  * tt/2  - aT> 


(A18) 


The  b.'s  of  Table  1 become 
3 

b1  = sec  0T;  b2  = b3  = b4  = b5  = 0;  bg  = 1 (A19) 

and 

■ • 

“l  = "ia  = “is  = tan  0T  + V’V  CA20) 

Equation  (A20)  for  0T  = 0 is  equivalent  to  Equation  (4.3)  on  page  150 
of  Reference  8. 
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